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Polarization-induced phase separation and re-entrant transition of two 
component-fermions in a one-dimensional lattice 
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By investigating the compressibility of one-dimensional lattice fermions at various filling factors, 
we study the phase separation and re-entrant transition within the framework of the Bethe ansatz 
method. We model the system using the repulsive Hubbard model and calculate compressibility as 
a function of polarization for arbitrary values of chemical potential, temperature, and interaction 
strength. For filling factors 0 < n < 1, compressibility is a non-monotonic function of polarization at 
all thermodynamic parameters. The compressibility reveals a phase transition into a phase-separated 
state for both low and intermediate temperatures at intermediate interactions as one increases the 
polarization. For certain filling factors, we Hnd the re-entrant transition into the mixed phase at a 
higher polarization. 


I. INTRODUCTION 


Ultracold fermionic atoms in optical lattices are cur¬ 
rently attracting a great deal of interest due to the 
possibility of impressive experimental simulation of rich 
physics associated with the strongly correlated condensed 
matter systems fli- While Feshbach resonance and 
laser intensities provide unprecedented control of atom- 
atom interactions in the optical lattice, laser interference 
phenomena provides the control of dimensionality. Two 
hyperfine states of Fermi atoms play the role of the up 
and down spins of the electrons. Unlike condensed mat¬ 
ter systems, the population of spin up and down parti¬ 
cles can be indraendently controlled by using a radio¬ 
frequency field |3l-H^. As a result, polarization can be 
maintained at any desired value between 0 and 100%. 
The term polarization here refers to a density imbalance 
of two hyperfine states. 


This study was motivated by the possibility of using 
ultra-cold atoms to engineer condensed matter systems. 
Condensed matter systems, such as transition metal ox¬ 
ides and rare-earth materials show collective and interac¬ 
tion dominated phenomena due to the electron-electron 
correlation effects. These phenomena such as the Mott- 
insulator transition and magnetism in strongly correlated 
materials is believed to be explained by the Hubbard 
model. The repulsive Hubbard model is the simplest 
model capable of explaining both metallic and insulating 
like behaviors, as well as the magnetic properties caused 
by electron correlation. The insulating states can be ei¬ 
ther a band insulator, caused by the Pauli exclusion of 
fermions, or a Mott-insulator, caused by strong on-site 
interactions. In the strong interaction limit, the local¬ 
ized magnetic properties depend on various parameters 
such as filling factors, orbital occupations, crystal field 
effects, and Hund’s coupling strength. 


Over the years, the Hubbard model has been the cen¬ 
ter of intense research as it captures the behavior of par¬ 
ent superconducting compounds and other magnetic ma¬ 
terials. The studies of phase separation for the Hub¬ 


bard model intensified after the experimental indication 
of phase separation of hole-rich and hole-poor regions in 
cuprate superconducting materials [l^. Strongly corre¬ 
lated electrons and holes are expected to play a key role 
in these materials and their phase separation is believed 
to hinder the superconductivity. If the Hubbard model 
is the correct model for the parent compound of super¬ 
conductors, can it be used to explain phase separation? 
This is the question that inspired us to study the phase 
separation of an exactly solvable one-dimensional model 
relevant to a flexible cold-atom experiment. It has been 
theoretically shown that there is no phase separation for 
two-dimensional bipartite lattices at any filling factors at 
finite temperatures 1^. In contrast, the phase separa¬ 
tion for the one dimensional Hubbard model is confirmed 
only close to a half filling in the presence of a critical 
magnetic field 1^. Finite-temperature phase separation 
for the one dimensional Hubbard model away from half 
filling has not been intensively investigated except for 
special cases Bin. 


In this paper we investigate the phase separation of 
one dimensional lattice fermions by Bethe ansatz (TBA) 
numerical method. We use the one dimensional Hubbard 
model as an effective model to describe the population- 
imbalanced two-hyperfine mixture in the optical lattice. 
We study phase separation by calculating the compress¬ 
ibility for various parameter regimes. By investigating 
the compressibility, we find phase separation at finite 
temperatures and at intermediate interactions as one in¬ 
creases the polarization. For some filling factors, we find 
a re-entrant transition into a mixed phase at a higher 
polarization. 


This paper is organized as follows. In section H, we 
discuss the geometry of the system and its connection 
to the one dimensional Hubbard model. In section HI, 
we briefly discuss our finite-temperature TBA calculation 
scheme. In section IV, we discuss the compressibility cal¬ 
culations and their connections to the phase separation 
and re-entrant transition into a mixed phase. We devote 
section V to discussion of the experimental connections 
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and we provide an experimental scheme to detect the 
phase separation. Finally in section VI, we draw our 
conclusions. 

II. THE MODEL: A ONE DIMENSIONAL 
OPTICAL LATTICE AND THE HUBBARD 
MODEL 

In general, a one dimensional optical lattice refers to 
an optical lattice generated by one set of laser stand¬ 
ing waves. The result of combined trapping and a pe¬ 
riodic potential gives a pancake-like shape of the sur¬ 
faces of constant potential. However, the geometry we 
consider here is generated by a three-dimensional opti¬ 
cal lattice where reduced dimensionality is achieved by 
freezing the atomic motion in the transverse direction. 
This can be done by operating two standing waves out 
of three mutually perpendicular laser standing waves at 
higher beam intensities. The higher intensities suppress 
tunneling in the transverse direction and create an ar¬ 
ray of one-dimensional lattice tubes. The dynamics of 
the atoms in one lattice tube can be modeled by the one 
dimensional Hubbard model given by, 

H = -t cj^Cja + UY (1) 

i(T icr 

The first term is the kinetic energy and is proportional 
to the tunneling amplitude t between lattice sites i and 
j = i + 1. The operator cj^(cicr) creates(destroys) a Fermi 
atom with hyperfine state denoted by pseudo-spin <j =t 
, 4 , (±1) at lattice site i. The second term describes the 
on-site interaction energy U. The density operator or the 
occupation number operator is riia = c\^Cia. Notice that 
< ij > indicates only the nearest neighbor pair of sites 
and we neglect tunneling beyond the nearest neighbors. 
The average chemical potential /x = (/x-j- -|- and the 
chemical potential difference h = (fj,^ — where /io- 

is the chemical potential of hyperfine state a. Here we 
neglect the confinement harmonic trapping potential and 
consider the lattice tubes are homogenous in space. The 
effect of trapping potential is discussed in section V. 


The tunneling amplitude and on-site interaction are re¬ 
lated to the complete set of Wannier functions Wn,iir) = 
ria localized at position fi with band index 

n, where a = x,y, z are the components of Cartesian 
coordinates. As the band gap becomes larger than U 
and temperature T, only the lowest band n = 0 is pop¬ 
ulated. For deep lattices, the lattice potential at site i 
can be approximated as a three-dimensional harmonic 
potential with frequency = 2Eii^/sZ/Ti, where SaEu 
is the laser intensity of the standing wave in the a di¬ 
rection. The recoil energy En = {fik)'^/2m is the ki¬ 
netic energy of an atom with mass m and the momen¬ 
tum hk of a single lattice photon. For deep lattices, 
taking WQ^a — a^) as a ground-state harmonic oscillator 
function with frequency Wq, the tunneling amplitude in 
the one dimensional geometry becomes t = J dxwQ(x — 
Xi)[—+ Vo{x)]wo{x — Xj). This is obtained from 
the Mathieu equation as t = A/y/nEiiSx^'^ {—2.^/3^), 
where Vo{x) = s^Ensiia^^kx) is the periodic potential 
generated by counter propagating lasers in the x direc¬ 
tion. The lattice constant d = A/2 is related to the laser 
wave length A, hence the wave vector k = 27r/A. The on¬ 
site interaction U = f dx\wQ{x)\^/m oc agy/sZ- 

Notice that the on-site interaction U can be repulsive or 
attractive depending on the free-space s-wave scattering 
length Os- In the present work, we consider a tight one 
dimensional geometry in the x direction with a positive 
U modeled by Eq. 1. Notice that the ratio t/U can easily 
be controlled by the laser intensity / oc Sa, of the counter 
propagating lasers in the x-direction. In our model, the 
laser intensities in the transverse directions that are pro¬ 
portional to Sy and Sz are maintained at higher intensi¬ 
ties so that the tunneling in the transverse direction is 
neglected. 


HI. THERMODYNAMIC BETHE ANSATZ 
METHOD 

Lieb and Wu have shown that the model presented in 
the previous section is exactly solvable in one dimension 
using the thermod yna mic Bethe-ansatz method [l^ . Fol¬ 
lowing Takahashi [lOl |20j , the thermodynamic potential 
per site is given by 


J 


n = eo - y - ksT^^ J po(fc) ln[I + -b J (To(A) ln[l-b ? 7 i(A)]dA|. 


( 2 ) 


1 f°° 

Po{k) = - -bcosA: / ai(A — sinA:)cro(A)dA (3) 

Stt j 


The energy per site here is given as eo = 2tl, and two 
distribution functions of fc’s and A’s are given by. 


— 00 
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1 /■°° 

cro(A) = — s{A — smk)dk. 

27r 7-00 

The two additional expressions introduced in the equa¬ 
tions are ai(x) = + IGa::^)] and s{x) = 

csc{2-kx/u)/u with u = U/t. The quantity I is related to 
the mth order Bessel functions Jm{x) through, 


The particle-hole ratios of k excitations and A excita¬ 
tions, ^(fc) and r]i{k) are obtained by an infinite set of 
nonlinear integral equations: 


/ = 



Jo{x)Ji{x) 
x(l + 


( 4 ) 

I 


(5) 

ln?7i(A) = s* ln[l-I-772(A)] — j s(A — sinA:) ln[l-|-^“ 1 (A:)] cosfcd/c (6) 

J —TT 

In 77](A) = s* ln[l-I-772(A)] — f s(A — sin A:) ln[l-I-^(fc)] cos/cdfc 

J —TT 


and for j >2, 


ln?7j(A) = s*ln{[l-H77j_i(A)][l-H77j+i(A)]} (7) 

In 77' (A) = s*ln{[l-H77'_i(A)][l-H77'+i(A)]}. (8) 


Here we use two integral functions given by s*/(A) = 
f^s(A — A')f{A')dA' and ko(^) = —2tcosk — 

At J^^dAs{A — sin A:) x Re^Jl — (A — ui/AY- The av¬ 
erage chemical potential and the chemical potential dif¬ 
ference are entered in the formalism through the grand 
potential fl, 


lim 

n—>CXD 


ln77„(A) 2h 
n ~ T’ 


(9) 


and 


lim = U-2^ ^ 

n->-oo n T 


( 10 ) 


In order to calculate the thermodynamic potential nu¬ 
merically, one has to cut off the set of infinite equations 
at a finite number j. We achieve this by following the 
numerical procedure proposed by Takahashi et al. [21 1. 
The infinite set of equations is truncated by replacing 
s(A) by 5{A)/2 at j > ric- Then the integral equations 
are converted into a set of matrix equations in which 
277c + 1 unknown functions are represented in terms of 
discrete points of k and A. These non-linear matrix equa¬ 
tions are then solved iteratively using Newton’s method 
for a given temperature (T), average chemical potential 


(/i), and chemical potential difference (h). The details 
of the numerical procedure can be found in Refs. HH- 
[^ . From the numerical solutions of the non-linear in¬ 
tegral equations, we first calculate the thermodynamic 
potential H using Eq. (2), and then the particle den¬ 
sity 77, = Ti-i- -|- rij, = —dn/dfi and the magnetization 
(the density difference of two hyperfine states, 77| — 77 ^) 
777 = did/dh follow. The compressibility is then calcu¬ 
lated numerically at a constant polarization P = m/n 
using the second derivative of the thermodynamic poten¬ 
tial with respect to the chemical potential [U, [1^ . 


IV. THE RESULTS: IDENTIFYING PHASE 
SEPARATION AND RE-ENTRANT 
TRANSITION THROUGH COMPRESSIBILITY 

We examine the stability of the mixed phase through 
the sign of compressibility. Negative compressibility indi¬ 
cates an instability of the mixed phase, where the system 
enters into a phase-separated state. 

Figure [ 1 ] shows the compressibility at a constant interac¬ 
tion strength and a constant finite temperature for vari¬ 
ous values of density. The compressibility is always pos¬ 
itive close to the densities of half filling and zero filling. 
However, away from these two limits, the compressibil- 
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FIG. 1: (color online) Compressibility of a one 
dimensional lattice fermion system at interaction 
strength U = 2t and temperature /3 = 30/t. As shown 
in inset, the density is approximately fixed. The 
negative compressibility indicates the instability of the 
mixed phase against the phase-separated state. Notice 
that the plotted quantities are dimensionless. 


ity becomes negative and then gets positive again as one 
increases the polarization. Notice the compressibility for 
filling factors n ~ 0.5 and n ~ 0.75 in the figure. The 
negative compressibility indicates the instability of the 
mixed phase, meaning that the system is phase separated 
into two different distinct phases corresponding to their 
pseudo-spins. A further increase of polarization causes 
the system to make a re-entrant transition into the mixed 
phase. In the mixed phase, both spin components coexist 
in the same region of space. This positive compressibility 
at higher polarization itself does not guarantee the sta¬ 
bility of the mixed phase over the phase separated phase. 
One has to compare the energies of phase-separated state 
and the mixed phase to determine the stability. The 
zero-temperature stability of the mixed phase at higher 
polarizations is justified in Ref. [3. This justification 
has been confirmed by comparing the ground-state ener¬ 
gies in both mixed and phase separated states using both 
weak and strong coupling approaches. We believe this is 
true even for finite temperatures. The comparison of fi¬ 
nite temperature energies of the phase-separated state 
and the mixed state is not trivial and these calculations 
are beyond the scope of the present paper. 

The zero-temperature instability of the mixed phase at 
higher polarizations has already been established within 
the bosonization theoretical frame work 0 . Bosoniza¬ 
tion theory suggests that phase separation occurs for 
U/t > 47r{sin[7r(n-|-m)/2] sin[7r(n —to)/ 2]}^/^. As shown 
in Fig [2 the mixed phase is stable only at higher den¬ 
sities, low polarizations, and low interaction strengths. 
The phase-separated state is stable at higher polariza¬ 
tions; however, unlike finite temperatures, the system 
does not make a re-entrant transition into the mixed 


0 



FIG. 2: (color online) Zero temperature phase diagram 
of one dimensional lattice fermions in polarization (P), 
on-site interaction ([/), and density (n) parameter 
space. The phase diagram is constructed from 
bosonization Theory. While the shaded region 
represents the phase separated state, the unshaded 
region represents the mixed state. 



FIG. 3: (color online) Compressibility of a one 
dimensional lattice fermion system at a constant 
temperature /3 = 30/t for various interaction strengths. 
See FIG. [T] caption for details. 


phase at zero temperature. 

The compressibility at various interaction strengths 
and temperatures is shown in Fig. |3] and Fig. |4] respec¬ 
tively. Here the compressibility is calculated at a de¬ 
sired polarization P by varying the chemical potential 
difference h and keeping the average chemical potential 
at a representative fixed value, /i = 3t. Notice that the 
mixed phase is stable for the entire range of polariza¬ 
tion at smaller and larger interactions. This can be jus¬ 
tified by the compressibility at the infinite interaction 
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FIG. 4: (color online) Compressibility of a one 
dimensional lattice fermion system at a constant 
interaction strength U = 2t for different temperatures. 

See FIG. [T] caption for details. 

and non-interacting limits. At the infinite interaction 
limit, the TEA equations can be solved to get analyti¬ 
cal results [ 2 ^. In the limit 17 —>■ 00, the polarization 
P = tanh(/3/i)/2 and the compressibility can be calcu¬ 
lated as K = 2/3cosh(/3/i)/4/J^/(fc)dfc, where f(k) = 
exp[/3(/r -I- 2tcosfc)]/(l -|- cosh(/3/i) exp[/3(^ -I- 2tcosfc)])^. 
As f{k) > 0 for all k values, the compressibility at the 
infinite interaction limit is always positive. This is intu¬ 
itive as the system can be considered as spinless fermions 
in this limit. On the other hand, in the limit 17 —>■ 0, no 
phase separation occurs as the system consists of non in¬ 
teracting fermions. Again, the negative compressibility 
at intermediate interactions suggests the phase separa¬ 
tion into two pseudo spin states. 

Consider the temperature dependence shown in FIG. SI 
The mixed phase makes a transition into a phase sepa¬ 
rated state and then makes a re-entrant transition into 
mixed phase at higher polarization for low temperatures. 
In contrast, the mixed phase is stable for higher tempera¬ 
tures (small f3) over the entire range of polarization. The 
high temperature expansion of the thermodynamic po¬ 
tential for the one dimensional Hubbard model has been 
carried up to the fourth order in /3 by Charret et al 12611 
and up to the sixth order in /3 by Takahashi et al [2l|. 
By using the 6th order expansion, we confirm the pos¬ 
itive compressibility at higher temperatures by an ana¬ 
lytic calculation. The high temperature expansion of the 
compressibility and the polarization up to the sixth order 
is given in the Appendix. 

Notice that compressibility is a non-monotonic func¬ 
tion of polarization for all temperatures and interactions. 
In contrast, compressibility is a non-monotonic function 
of the interaction parameter only for larger polarizations. 
However, as evident from the Fig. [3l compressibility is a 
monotonic function of temperature for the entire range 
of polarizations. 

It is worth mentioning that a small density imbal¬ 



FIG. 5: Spatial variations of the atom density 
n{z) = n^{z) + n^^z) [black curve in panel (a)], 
magnetization m(z) = n^{z) — n 4 ,(z) [gray curve in 
panel (a)]) , and polarization 
P{z) = [n^{z) - ni{z)]/[n^{z) + ni{z)] [panel (b)]. We 
define the scaled length z = z^mixP 12, where a; is the 
one dimensional trapping frequency. We fixed the 
on-site interaction (C7 = 2t) and the inverse temperature 

ance can be induced in condensed matter electronic sys¬ 
tems by applying an external magnetic field. Thermo¬ 
dynamic properties of such one dimensional systems are 
thoroughly discussed in Ref. . Though finite temper¬ 
ature compressibility as a function of polarization is not 
discussed in there, special attention has been given to 
the ground-state properties such as susceptibility, mag¬ 
netization, and densities ( 2 ^ . 

V. CONNECTIONS TO EXPERIMENTS 

Recent progress in experimental techniques with ultra¬ 
cold atoms, such as single-site detectio n . noise 

correlations Ej, Bragg scattering [32|. and in situ 
imaging in the lattice scaling [ 2 ^, allows one to probe 
the density variations in cold-atom experiments. For the 
case of equal-population two-component fermions on a 
three-dimensional cubic lattice, the compressibility has 
already been measured . 

Though we neglected it in this study, the underlying 
harmonic trapping potential present in all cold gas ex¬ 
periments causes the density to vary across the lattice. 
See FIG. [3 By combining the TBA solutions with the 
local density approximation (LDA), we then extract the 
local density n(z), magnetization m(z), and then polar¬ 
ization P{z). In LDA, the external trapping potential 
Vi = muP'z^ 12 at site i is related to the local chemical 
potential through the relation fj,i = fig — Vi, where w is 
the one-dimensional trapping potential, /io is the central 
chemical potential and z = id, with lattice constant d the 
spatial coordinate. As shown in Fig.[3 the density mono- 
tonically decreases, while polarization monotonically in¬ 
creases from the center to the edge of the trap [i^. This 
trapping potential induced inhomogeneity allows both 
mixed-phase and phase-separated states to exist simul¬ 
taneously inside the trap. At the center of the trap, the 
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density is higher and the polarization is lower. On the 
other hand, the polarization is higher and the density is 
lower at the edge of the trap. As a result, the mixed- 
phase should exist at the center and at the edge of the 
trap. However, depending on the density, the phase sep¬ 
arated state can exists in the middle (not the center) of 
the trap. Therefore, by adjusting the total density in 
the trap, any polarization induced phase separation and 
re-entrant transition can be investigated experimentally 
with currently available experimental techniques. 


We find that compressibility is a non-monotonic function 
of polarization. At filling factor 0 < n < 1, with low tem¬ 
peratures and intermediate interactions, compressibility 
becomes negative, indicating instability of the mixed- 
phase state towards the phase separated state. For some 
parameters at higher polarizations, compressibility be¬ 
comes positive again indicating a re-entrant transition in 
to a mixed phase. These phase-separation and re-entrant 
transitions can be detected by using currently available 
experimental techniques. 


VII. CONCLUSIONS 


VIII. APPENDIX: HIGH TEMPERATURE 
EXPANSION OF THE COMPRESSIBILITY 


In conclusion, we considered two-component Fermi 
atoms in a highly tunable optical lattice to study the 
phase separation of fermions in one dimension. We have 
calculated the compressibility of one dimensional lattice 
fermions using the thermodynamic Bethe ansatz method. 


In this appendix, we provide compressibility and po¬ 
larization up to the sixth order in j3. Using the high- 
temperature expansion of the thermodynamic poten¬ 
tial 21, 2^, we find the compressibility. 



C //32 


32384^ 


^5 

1536 




Kq 


30720 


/3® + 0(/9"), 


where the coefhcients at each order are given by 


( 11 ) 


K3 = -4h2-8t2-(17-2/4)2, 


( 12 ) 


ka = t/(48t2 -/ 13172 - 4817/4 -h 48/42), 


(13) 


K5 = 32/i^ -h 192t4 -h 96^2(17 - 2/4)2 _ _ 2^f{7U'^ + 817/4 - 8^^) + 4h^[15U^ + 2t^{A8 + U^) - 4817/4 -h 48/42],(14) 

and 

K6 = -U[-2A0h‘^ + 2880t'^ + 107U^ + 860h^(U - 2^if - 7601/^ + 2l2W‘^p? - 2720Un^ + 1360/4^ (15) 

-h240t2(lH72 - 4017/4 -h 40/42)]. 


The sixth-order expansion of polarization becomes. 


hp h{U - 

'2 4 





46080 


/3® + 0(/3"), 


where the coefficients of higher orders are. 


(16) 


^3 = -hijr^ + 6 ^ 2 ), 


(17) 


Pa = -h[Ah^ + 12^2 -h (17 - ^if]{U - fi) 


(18) 


P 5 = /4{16l4^ -f 5t2[96t2 -/ (12 -h U^){U - 2/4)2] -h 40h^[8t^ - 3(17 - nf)]}, 


(19) 
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and 


Pe = h{ilQh^{U - Ai) + 20(/i^)[16(P - nf + 3^2(160^ +U^ - 160/^ - 2PV)] + 3[2880t^(P - m) + 32(P - (20) 

+5t2(p5 - 6PV + 576 Pm^ - 256/^3 - 8PV(60 + ^Ji^) + 4p3(38 + 3//^))]}. 
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